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ABSTRACT 

We present a re-appraisal of the temperatures of Red Supergiants (RSGs) using their optical and 
near- infrared spectral energy distributions (SEDs). We have obtained data of a sample of RSGs in 
the Magellanic Clouds using VLT-^XSHOOTER, and we fit MARCS model atmospheres to different 
regions of the spectra, deriving effective temperatures for each star from (a) the TiO bands, (b) line- free 
continuum regions of the spectral energy distributions (SEDs), and (c) the integrated fluxes. We show 
that the temperatures derived from fits to the TiO bands are systematically lower than the other two 
methods by several hundred Kelvin. The TiO fits also dramatically over-predict the flux in the near- 
IR, and imply extinctions which are anomalously low compared to neighbouring stars. In contrast, the 
SED temperatures provide good fits to the fluxes at all wavelengths other than the TiO bands, are in 
agreement with the temperatures from the flux integration method, and imply extinctions consistent 
with nearby stars. After considering a number of ways to reconcile this discrepancy, we conclude that 
3-D effects (i.e. granulation) are the most likely cause, as they affect the temperature structure in the 
upper layers where the TiO lines form. The continuum, however, which forms at much deeper layers, 
is apparently more robust to such effects. We therefore conclude that RSG temperatures are much 
warmer than previously thought. We discuss the implications of this result for stellar evolution and 
supernova progenitors, and provide relations to determine the bolometric luminosities of RSGs from 
single-band photometry. 

Subject headings: Stars: fundamental parameters. Stars: atmospheres, Stars: late-type, Stars: mas- 
sive, supergiants. Stars; evolution 



1. INTRODUCTION 

Red Supergiants (RSGs) are a post main-sequence 
(MS) phase of s tar s wit h masses ^8-30M f7) (e.g . 
Mevnet fc M aedei^ 120001 : lEldridge et all 120081 : 
Brott et al. 2011) . Their lum inosities (^ lO^-^-^-^L© 
Humphrevs fc DavidsonI Il979f ) rival those of globular 
clusters and dwarf galaxies, and dominate the light 
output of their host galaxies at near-infrared (IR) 
wavelengths. For most stars in this mass range, the 
RSG phase is thought to be the evolutionary stage 
immediately preceeding core-collapse supernova (SN), 
though searches of pre-explosion images are yet to find 
a RSG progen itor with an inferred initial mass gre ater 
than ~2OM0 (iSmartt et al.ll2009l : HVaser et al.ll2012[ ). 

In determining the bolometric luminosities of RSGs, 
knowledge of the temperature scale is crucial, since 
for the typical temperature range of RSGs (Toft'^3400- 
4500K) the optical bolometric correction can vary by al- 
most 2 magnitudes (jLevesaue et al.|[2005f) . Therefore, in 
order to accurately convert the flux in a given photo- 
metric pass-band to the total flux, we must know the 
relationship between the star's observed properties (e.g. 
photometric colours, spectral type) and its Toff. This 



problem is complicated further by the fact that extinc- 
tion Ay ^ whether interstellar or circumstellar, can make 
the star appear cooler, introducing a large degree of de- 
generacy between Ay and Toff. 

To measure the effective temperatures of RSGs, a num- 
ber of different strategies have been employed, though all 
methods ultimately require a knowledge of the radiative 
transfer through the star's atmosphere. One method is 
to measure the angular diameters of stars through inter- 
ferornetry or Lunar occul t ation s (iLedllQTOl: iDvck et al.l 
119961 : iMozurkewich et al.l l2003t Ivan Belle et al.l l2009| l. 
Since the definition of Tefi comes from L^, = AirRlaT^Q, 
where L^, and i?^ are the stellar luminosity and radius 
respectively, by measuring the apparent sizes and fluxes 
one can determine Toff. 

This method does have a number of problems how- 
ever. Firstly, the measured size of the star is very sensi- 
tive to the observed wavelength, with the star appearing 
larger in spectral regions of high opacity. For example, 
interferometric observations in the near-IR have shown 
that the apparent size of the star may increase by up 
to 50% around the wavelengths of str ong molecular lines 
compared with the nearby continuum (jPerrin eraI1l200l 
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Ohnaka et all [200l Pfsu^ [2OOOI: lOhnaka et all [2011 
Chiavassa et alJ 12010( 1. Therefore, in order to interpret 
the measurements of angular diameters, model stellar at- 
mosphere are required. Early work typically used black- 
bodies, which are well-known to be a poor representation 
of RSGs, particularly in the optical. Later work has typ- 
ically used plane-parallel Kurucz models, again known to 
be inadequate for RSGs. 

Secondly, since one needs to know the flux of the star, 
one also needs to know the foreground extinction. Typ- 
ically it is assumed that the target star has the same 
foreground extinction as the other stars in its host clus- 
ter or association. However, this is likely to be an under- 
estimate, since RSGs can produce their own circumstel- 
lar dust through their winds, provi ding up to a magni- 
tude of extra visual extinction (e.g. iSchuster et al.|[2006l : 
Ide Wit et al.l[20 08). A combination of these two effects 
could explain why near-IR angular diameter measure- 
ments provide warmer effective temperat ures (i.e. smaller 
stellar radii) than those in the optical (jvan Belle et al.l 
[20091 

Other investigations into the temperatures of RSGs 
have loo ked at the i r spe ctral energy distributions 
(SEDs). iDvck et al.l (|1974D defined a near-IR 'colour 
temperature ', work which was built upon by iFloweij 
([T975I [T97I . This work was however calibrated this 
against angular diameter measurements taken at 7100A, 
which is at the centre of a deep absorption feature, and 
so may have over-estimated the radii. Later work fit 
data with s ynthesiscd spect ra from pla ne-parallel models 
([White fc Win g 1978; Scar gle fc Strecker 197 9). though 
these models lacked detailed molecular opacities. 

The most contemporary m easurement of the RSG tem- 
perature scale is presented in lLevesque et al.l ([20051 [20061 
hereafter LOS and L06). These authors obtained spec- 
trophotometry of samples of RSGs in the Galaxy and 
the Magellanic Clouds in the B VRI region of the spec- 
trum. They then fitted this region of the spectrum, which 
is dominated by the TiO absorption bands which de- 
termine the spectral classifica tion of these stars, wit h 
MARCS model atmospheres ([Gustafsson et al.l l2008f ). 
which are 1-D spherical hydrostatic models, important 
for stars with extended atmospheres. The temperature 
scale they found for Galactic stars was slightly warmer 
than was previously thought. These warmer tempera- 
tures brought Galactic RSGs into closer agreement with 
the predicted Haya.shi li mit of Geneva stellar models^ 
([Mevnet fc Maedeil[2000[ ). 

There are, however, known issues and anomalies with 
the LOS Teff scale as well. These authors found that 
the temperatures derived from fitting the optical spectra 
were systematically offset from those expected from the 
V — K colours, indicating that their model fits were not 
reproducing the correct near-IR flux. They also found 
a number of very cool stars in the SMC, whose tem- 
peratures were so low that they were apparently in vi- 
olation of hydrosta tic eauilibriun i (iLevesaue et aP 120071 : 
iMassev et al]|2007l ). In addition. iLancon et al.l ([20071 ) fit 

^ It should be noted that the location of the Hayashi hmit in 
models of massive stars is governed by the convective mixing length 
parameter, whieh is typically assumed to be the same as for the 
Sun. Observations of globular clusters have shown that a uniform 
mixin g length may not b e appropriate for all stars, especially cool 
stars jFerraro et al.1120061 '). 



the optical and near-IR spectra of a sample of Galactic 
RSGs with PHOENIX models, finding problems simulta- 
neously fitting the strengths of various molecular bands, 
particularly for the most luminous stars. 

The reasons for these discrepancies could lie in the de- 
tailed radiative transfer that governs the strengths of the 
TiO bands, which are used as temperature indicators by 
LOS and L06. As will be discussed later, the strengths of 
these bands are sensitive to a number of factors, such as 
the temperature structure of the atmosphere and metal- 
licity. Such factors may make the TiO bands unreliable 
temperature indicators. 

Here we explore an alternative method of deriving the 
temperatures of RSGs. As in LOS and L06, we fit spec- 
tra with MARCS model atmospheres. However, in ad- 
dition to looking at the strengths of the TiO bands, we 
also investigate the implied by fits to the optical- 
infrared SED using regions which are free of the deep 
molecular absorption seen in the optical regions. We do 
this on a sample of stars taken from the two Magellanic 
Clouds, which should have relatively low foreground ex- 
tinction, and where the discrepancies in the Levesque 
et al. temperature sc ale appear to be most conspicuous 
([Levesque et al.ll2007l ). 

We begin in Sect. [2l with a description of our observa- 
tions and data reduction. In Sect. [3|we discuss our anal- 
ysis techniques, and present our results in Sect. [4| which 
indicate a systematic error in the LOG temperature scale. 
Reasons for this discrepancy, and possible remedies, are 
discussed in Sect. [Sj as well as the implications of our 
results. We conclude in Sect. [S] 

2. OBSERVATIONS & DATA REDUCTION 

We have obtained observations of several stars 
in the LMC and S MC using VLT+XSHOOTER 
(iD'Odorico et al.l[2006[) under ESO programme number 
088.B-0014(A) (PI B. Davies). The sample of stars, 8 
in the LMC and 10 in the SMC, were selected from L06. 
We aimed to sample the full distribution of spectral types 
observed in each galaxy. We also deliberately excluded 
stars which may be 'extreme' objects, such as analogues 
of the Galactic RSG VY CMa. We did so by avoiding 
stars which had luminosities in excess of 10^'^Lq in L06. 

The stars were observed in nodding ABBA mode with 
four exposures per star, and a randomized jitter at each 
position on the slit. For each of the instrument arms - 
UVB, VIS, and NIR - we used the 5.0"sfits to minimize 
slit losses and obtain accurate spectrophotometry. The 
spectral resolution was therefore set by the seeing, which 
for seeing of 1" was roughly R = A/ A A ~ SOOO at all 
wavelengths. The precise value of R was determined at 
the analysis stage (see Sect. [3|). Integration times are 
listed in Table [1] and were chosen to achieve a signal-to- 
noise ratio (SNR) of at least 50 per resolution element at 
all wavelengths greater than 400nm. Since the UVB and 
VIS arms take much longer to read-out than the NIR 
arm, and each arm operates independently, extra time 
could be spent integrating in the NIR. 

Flux standard stars were observed each night and to 
correct for the atmospheric absorption in the NIR, tel- 
luric standard stars of spectral type late-B were observed 
within one hour of each science target. In general observ- 
ing conditions were good on each night, and the seeing 
was always <1". The standard suite of XSHOOTER cal- 
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TABLE 1 

Observation data for the stars in our sample. The errors on the synthetic photometry are ±0.1mag for the optical data 

AND ±0.04MAG for, the INFRARED DATA. 



Star 






RA DEC 
(J2000) 






Obs date 


£UVB 
(sec) 


^VIS 
(sec) 


^NIR 
(sec) 


B 


V 




I 


J 


H 


Ks 


SMC 


011709 





48 


46.32 


-73 


28 


20, 


.7 


2011- 


-10-13 


60 


10 


50 


14.43 


12.60 


11, 


,63 


10 


,76 


9.63 


8.84 


8.61 


SMC 


013740 





49 


30.34 


-73 


26 


49, 


.9 


2011- 


■12-06 


80 


12 


70 


15.56 


13.78 


12, 


.76 


11 


.86 


10.61 


9.75 


9.49 


SMC 


020133 





51 


29.68 


-73 


10 


44, 


.3 


2011- 


-12-06 


60 


10 


50 


14.84 


12.86 


11, 


.76 


10 


.73 


9.39 


8.56 


8.27 


SMC 


021362 





51 


50.25 


-72 


05 


57, 


.2 


2011- 


■12-05 


80 


12 


70 


14.91 


13.02 


12, 


,00 


11, 


.09 


9.91 


9.05 


8.80 


SMC 


030616 





54 


35.90 


-72 


34 


14, 


.3 


2011- 


■12-06 


60 


10 


50 


14.50 


12.67 


11, 


.65 


10, 


.69 


9.48 


8.68 


8.45 


SMC 


034158 





55 


36.58 


-72 


36 


23, 


.6 


2011- 


■12-06 


80 


12 


70 


14.87 


13.01 


12, 


.01 


11, 


.14 


10.04 


9.26 


9.03 


SMC 


035445 





55 


58.84 


-73 


20 


41, 


.4 


2011- 


■12-06 


100 


15 


100 


14.59 


12.91 


12, 


,00 


11, 


.20 


10.09 


9.31 


9.09 


SMC 


049478 


1 


00 


41.56 


-72 


10 


37, 


.0 


2011- 


■12-03 


50 


5 


30 


14.23 


12.40 


11, 


,35 


10, 


,36 


9.20 


8.42 


8.15 


SMC 


050840 


1 


01 


15.99 


-72 


13 


10, 


.0 


2011- 


■12-06 


60 


10 


50 


14.62 


12.72 


11, 


.66 


10, 


.69 


9.48 


8.64 


8.38 


SMC 


057386 


1 


03 


47.35 


-72 


01 


16, 


.0 


2011- 


■12-30 


80 


12 


70 


14.22 


12.71 


11, 


,83 


11, 


,07 


10.01 


9.24 


9.01 


LMC 


064048 


5 


04 


41.79 


-70 


42 


37, 


.2 


2011- 


■12-05 


60 


10 


50 


15.09 


13.15 


11, 


,94 


10, 


,69 


9.42 


8.55 


8.19 


LMC 


067982 


5 


05 


56.61 


-70 


35 


24, 


.0 


2011- 


■12-05 


60 


10 


50 


14.80 


12.84 


11, 


.73 


10, 


.61 


9.27 


8.40 


8.11 


LMC 


116895 


5 


19 


53.34 


-69 


27 


33, 


.4 


2011- 


■12-06 


60 


10 


50 


14.49 


12.62 


11, 


.54 


10, 


.44 


9.17 


8.39 


8.08 


LMC 


131735 


5 


23 


34.09 


-69 


19 


07, 


.0 


2011- 


■11-29 


80 


12 


70 


14.33 


12.60 


11, 


,73 


10, 


.98 


9.87 


9.15 


8.95 


LMC 


137818 


5 


27 


14.33 


-69 


11 


10, 


.7 


2011- 


■11-29 


60 


10 


50 


15.34 


13.38 


12, 


.03 


10, 


.65 


9.34 


8.52 


8.18 


LMC 


142202 


5 


28 


45.59 


-68 


58 


02, 


.3 


2011- 


■12-06 


50 


5 


30 


14.25 


12.30 


11, 


.25 


10, 


.19 


8.78 


7.95 


7.64 


LMC 


143877 


5 


29 


21.10 


-68 


47 


31, 


.5 


2011- 


■10-13 


50 


5 


30 


14.04 


12.13 


11, 


,13 


10, 


,21 


9.00 


8.28 


8.01 


LMC 


158317 


5 


33 


44.60 


-67 


24 


16, 


.9 


2011- 


■10-13 


80 


12 


70 


15.12 


13.13 


12, 


,01 


10, 


,97 


9.64 


8.81 


8.48 



ibration frames used in the data reduction process were 
tak en at the beginning and end of the night (for details, 
see iGoldoni et aLl[20M) . 

The initial steps of the reduction process were 
done using the X SHOOTER data reduction pipeline 
(jModigliani et alJ 12010;) . These steps included subtrac- 
tion of bias and dark frames, flat-fielding, order extrac- 
tion and rectification, flux and wavelength calibration. 
The accuracy of the wavelength solution was checked by 
measuring the residuals of the wavelengths of the arc 
lines. The root-mean-square of these residuals was found 
to be below '^0.1 pixels (^Ikm s^^) for all targets. To 
check the flux calibration we applied the response func- 
tion to the flux standard and compared with the cali- 
bration spectrum. We found residuals of order 1%. The 
absolute accuracy of the flux calibration is discussed in 
Sect. 12.21 When stitching the spectra of the three differ- 
ent arms together, we clipped the UVB arm at 563nm 
(to avoid the flux dropout due to the first dichroic), the 
VIS arm between 558-lOlOnm, and the NIR arm below 
2350nm. The overlapping regions had consistent fluxes 
to within 1%, and the flux at these regions was calculated 
as the mean of the two overlapping spectra. 

The spectra of the science targets and the telluric 
standards were then extracted from final rectified two- 
dimensional orders. The NIR spectra of the telluric 
standards were found to have some non-gaussian noise, 
mainly from bad pixels and cosmic-ray hits. To identify 
and remove this noise a new spectral extraction algo- 
rithm was written. Briefiy, bad pixels were fiagged by 
comparing the spatial profile at a given spectral chan- 
nel with the median spatial profile of the neighbouring 
3 spectral channels either side. These pixels were then 
replaced with the corresponding value in the median pro- 
file. Typically, around five such pixels were found and 
replaced per spectral order. 

2.1. Telluric correction 

The method of removing the telluric absorption fea- 
tures depended on the observing arm. For the UVB arm, 
telluric features are only present at the blue end of the 
wavelength range, where the SNR of our science data was 
very low due to their late spectral types. For this reason, 
no telluric correction was performed for the UVB arm. 



In the VIS arm, telluric features are non- negligible, 
with prominent telluric features at 686nm, 717nm, 
759nm, 815nm, and 900nm. To correct the spectra in 
each of these regions, we used a synthetic telluric spec- 
trum, rather than the standard stars which have many 
intrinsic spectral features at these wavelengths, such as 
the Balmer hydrogen lines. The synthetic spectrum was 
optimized to best match the telluric absorption in an it- 
erative process. In each region the relative shift of the 
synthetic spectrum with respect to the science data was 
found by cross correlation. The strength of the absorp- 
tion features and effective resolution was then tuned to 
give the best cancellation. This whole process was then 
repeated. 

In the NIR arm, a similar process was employed, only 
this time the spectrum of the telluric standard star was 
used. As we wish to retain flux calibration of the science 
data, we flrst flattened the telluric spectrum to remove 
the intrinsic continuum slope, and removed the Brackett 
hydrogen lines by fitting voigt profiles. The stellar con- 
tinuum was corrected by fitting a black-body curve to the 
regions free of telluric absorption (the same regions used 
by the XSHOOTER pipeline in the fiux-calibration step 
of the reduction). The temperature of the model black- 
body was chosen to ensure that the ratio of its spectrum 
to the fiux-calibrated telluric standard star was closest 
to unity in all spectral regions tested. The accuracy of 
this fit was better than 5% in all tested regions of the 
spectrum. 

The remaining spectrum, which was a pure measure 
of the telluric absorption, was then tuned to match the 
science data. To do this, the data was first split into 
chunks of width lOnm. For each chunk, the optimal shift, 
absorption strength and resolution was found as in the 
VIS arm. We then assumed that the absorption strength 
would be uniform across the whole spectrum, while the 
shift and resolution would vary linearly as a function of 
wavelength. The absorption strength at all wavelengths 
was taken to be the median of all those measured in 
each of the spectral chunks. The optimized resolution 
and wavelength shift as a function of wavelength were 
obtained from linear fits to the results from the separate 
spectral chunks. 
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2.2. Synthetic photometry 

To check the accuracy of our photometric caHbration, 
and for comparison with future and past measurements, 
we derived synthetic photometry for all our stars. We 
convolved our XSHOOTER spectra with Bessell BVRI 
filter profiles, as well as with those of the 2MASS JHKg 
filters. 

The dominant sources of error in the synthetic photom- 
etry are the accuracy of the spectrophotometric response 
function, and for bands which have significant telluric ab- 
sorption, the quality of the telluric correction. Tests of 
the response function, calculated by the XSHOOTER re- 
duction pipeline, showed that it should be accurate to a 
few percent, providing that the observing conditions did 
not change between observing the science targets and 
the flux standards. Similarly, imperfections in the tel- 
luric correction should only account for a few percent of 
the total measured flux. 

To check the photometric precision and accuracy em- 
pirically, we compared our synthetic photometry with 
that of the M agellanic Clouds Photometric Survey 
(iZaritskv et al. 2002. 2004)and the 2MASS point-source 
catalogue ( Cutri et al.ll2003l) . Though many of the stars 



in our sample are known to be variable, the standard de- 
viation of the differentials between our photometry and 
that in the literature can provide an upper limit to the 
experimental uncertainties, whilst also giving an indica- 
tion of any systematics. We find that the optical BVI 
photometry has la standard deviations of ±0.2mags, 
while for the 2MASS JHKs photometry it is ±0.08mags. 
The systematic offsets are consistent with being zero. 
Given the variable natu re of RSGs. which can have am- 
plitudes ^0.5 mag (e.g. iKiss et al.|[2006t iLevesque et al.l 
l2007t lYang fc Jiang 201^ . a reasonable conservative es- 
timate of our photometric errors is around half these val- 
ues, or ibO.lmags for the optical data and ±0.04 for the 
infrared. 

3. ANALYSIS 

For this study we have computed a whole new grid of 
model atmospheres. The se atmospheres were ge nerated 
using the MARCS code (|Gustafsson et al.ll2008l) . which 
operates under the assumptions of LTE, spherical sym- 
metry and hydrostatic equilibrium. The grid of models 
is four-dimensional, computed with a range of metal- 
hcities (log(Z/Zo) = [Z]), gravities (logg), effective 
temperatures (Tcff) and microturbulent velocities (^). 
All model atmospheres were computed at ^=2 km s^^, 
with synthesized spectra computed at 2 km s~^ and 
5 km s~^ . The r elativ e abundance ratios were taken from 
iGrevesse et al.l ()2005() . Synthetic spectra were computed 
from the model atmosphere s using the turbospectrum 
code (jAlvarez fc Plezlll998l: lPlej|2012[) . with a spectral 
resolution of i? = 500, 000 between 250-2500nm. The 
chemical composition was scaled from Solar at [Z]=-1.5 
and -1-0.5 in steps of 0.25dex; Tcff between 3400 and 
4000K in steps of lOOK, with further models at 4200K 
and 4400K; log g between -0.5 and -1-1.0 in steps of 0.5 (in 
cgs units). All models were computed with an adopted 
stellar mass of M,t=15MQ. Though RSGs may have 
masses between ^8-25M0, the pressure scale height re- 
mains largely unchanged througho ut this mass range (see 
discussion m iDavies et aLll2010[ ). Within this grid of 



models, we logarithmically interpolated model spectra at 
metallicities appropriate for the LMC and S MC, specifi- 
cally [Z]=-0.65 and -0.4 ()Trundle et al.ll2007D . with mod- 
els ±0.25dex either side of these to investigate the effect 
of variations in metallicity. We also logarithmically in- 
terpolated the temperatures onto a finer grid of spacing 
20K. 

To derive the temperatures of the stars, we employ 
three different techniques, and compare the results of 
each. Below we describe the methodolgy of each. 

3.1. The TiO method 

Firstly, we extract a region of the spectrum around 
the TiO bands between 500-800nm. We compare the 
fiux in this spectral window with those of a subset of 
models. We initially selected models at the appropri- 
ate metallicity for the star's host galaxy (see previous 
section), log g=0.0, and ^=2km a~^. The choice of grav- 
ity is motivated by the typical masses, luminos ities and 
radii of RSGs (see e.g. lMevnet fc MaedeHl2000t) . Micro- 
turbulence^ parameters are typically between 2-4 km 
(jCunha et a l. 2007). In practice, varying metallicity, 
log (7 and f was found to have little effect on our results 
(see below). The parameters we are therefore allowing 
to vary are the temperature Tes and the extinction Ay, 
with the latter being computed from the reddening law of 
iGordon et al.l ([2003). Broadly, Tcff affects the strengths 
of the TiO absorption, while Ay affects the overall slope 
of the spectrum. We find the best fitting model by com- 
puting the summed over all spectral pixels for each 
value of Toff and Ay ■ We then fit for the two parameters 
simultaneously by finding the model with the minimum 
in Tcs-Av space. This method weights all TiO bands 
within the wavelength range equally. 

The errors in Tgg and Ay are measured from the max- 
imal values of those models which have values within 
0.5 of the best fitting model - this was found to be the 
point at which the quality of the fit becomes notice- 
ably poorer. We also explored the best fitting models 
with [Z]=±0.25dex, with logg=-0.5 — hl.Odex, and with 
^=5 km s~^. If these fits provided values of Tcff and Ay 
which were outside the errors of the initial fits, then the 
errors on these values were replaced by the best fits of 
these other subgrids. In general though, these differences 
were small (<50K). 

3.2. The SED method 

The second method of measuring the temperatures was 
to apply the same technique, but to the whole of the SED 
except those regions dominated by molecular absorption, 
since the strengths of these bands may be sensitive to 
the layer at which they form in the upper atmosphere. 
We avoid the whole of the BVR spectral region, the deep 
absorption bands of TiO and VO, the several narrow CO 
bands at ■~1.6/.tm and the band-heads >2.3/im, and the 
CN band at 1.1/im. 

Before fitting, spectra were first smoothed by a box- 
car filter of width 10 pixels in order to wash out narrow 
spectral features, as it is the broad shape of the SED that 
we are using as a diagnostic of the temperature. As with 
the TiO method, we used a x^-minimization technique to 
find the best-fitting v alues of Tcff and Ay We again used 
the reddening law of iGordon et al.l (|2003( ). and explored 



The temperatures of Red Supergiants 



5 



the effects of varying gravity, abundance and microtur- 
bulence in computing our error estimates. We also exper- 
i mented with dif f erent extinctio n laws, specifically th ose 
of lCardelli et"all (fToSOh and Rick e fc Lebolk^ psl . as 
well as values of Ry between 2-6 for the Cardelli et al. 
law. We found that varying the extinction law had only 
a very minor effect on our derived temperatures, result- 
ing in differences to our best-fit temperatures which were 
less than ±50K. 

3.3. The Flux Integration method 

The third method uses the definition of the effective 
temperature, 

aT^s = ^ / F,dX (f ) 

Jo 

by a wavelength integral over the stellar surface flux, 
where a is the Stefan-Boltzmann constant. This method 
takes advantage of XSHOOTER's wide wavelength range 
which allows to measure the effective temperature defin- 
ing flux integral directly provided that the reddening and 
angular diameter are known. We call this method the 
flux integration method (FIM). 

A straightforward application of the method is to as- 
sume a fixed value of interstellar reddening E{B — V). In 
this case we can determine angular diameter and Toff by 
a simple iterative procedure. W e deredden the obser ved 
SED using the reddening law of iGordon et al.l ()2003[ ). a 
reddening E(B — V), and an initial estimate of Toff 7 for 
instance the value obtained from the TiO fit. We then 
compare the observed flux with the model atmosphere 
flux at a particular wavelength (see below), and using 
the initial T^s estimate we then calculate the angular di- 
ameter. With this angular diameter we can then turn the 
observed flux at all wavelengths into the stellar surface 
flux, calculate the flux integral, and hence a new effec- 
tive temperature. The process is then iterated until the 
input and output temperatures converge, which usually 
occurs within four to five iterations. 

We choose to normalise the model and observed fluxes 
using the wavelength range 2.17-2.20/im, for the follow- 
ing reasons. Firstly, reddening at these wavelengths is 
minimal. Secondly, at these wavelengths we are close 
to the Rayleigh-Jeans domain of the SED. Here, one 
can analytically show that the process must converge as 
long as the temperature dependence of the model atmo- 
sphere flux at the wavelength Aq of the normalization 
pmodei ^ rpx^ goes with an exponent x smaller than 
four. Finally, the fact that the molecular line opacity 
at this wavelength range is weak reduces the uncertainty 
of the model atmosphere fluxes, and of the height in the 
atmo sphere at which the molecules form (jPerrin et al.l 
l2004f ). Moreover, at these wavelengths the line and con- 
tinuum form at roughly the same depth in the atmo- 
sphere as that which deflnes the photospheric radius in 
the MARCS model atmospheres. 

The advantage of this method is that is only very 
weakly model dependent through the iterative estimate 
of the angular diameter at the wavelength Aq and that it 
uses directly the deflnition of Toff- Finally, we note that 
this method is not new. A similar a pproach has been 
appli ed for hot stars (see for instance iRemie fc LamersI 
Il982f ) and also for cool stars (|BlackweU fc ShalhslllQTTD 



TABLE 2 

Results from the TiO and SED methods of estimating Tj.fi. 
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4. RESULTS 

4.1. TiO temperatures versus SED temperatures 

In Fig. ini we show examples of the fits from the TiO 
and SED techniques to two of the stars in the sample. 
Typically, the TiO fit is excellent in the region of the TiO 
bands (left panels). The fitted temperatures and extinc- 
tions are in most cases consiste nt to within ±100K o r 
itO.lmag with those obtained bv iLevesque et al.] (|2006( ). 
who used similar models with 'by-eye' fitting methods. 
Our errors are asymmetric due to the non-linear relation 
between model temperature and TiO strength, while for 
higher TiO temperatures the errors tend to be larger, 
due to the diagnostic TiO bands becoming weaker. The 
SMC stars have less TiO absorption for a given tem- 
perature due to the reduced metal abundance. This ef- 
fect explains the average spectral type of RSGs shift- 
ing to earlier subt ypes in lower metallicity environments 
(|Ehas et all 11981 . 

Longward of ~800nm, however, the fit is poor (right 
panels). Specifically, the model with the T^s of the 
TiO fit overpredicts the fiux at wavelengths greater than 
^900nm, as well as a rise in flux in the i/-band at 
~1.6/xm, which we call the ^H-hump\ It is caused by 
a minimum in the H" opacity, but is suppressed at high 
metallicities by molecular absorption, mostly CO. We do 
not see the H-hump feature in our XSHOOTER data. 
This discrepancy is not an artifact of the fltted extinc- 
tion: the grey lines in the right panels of Fig.[9]show mod- 
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Fig. 1. — Comparison of the temperatures derived from the TiO 
band fits and from the SED fits. 
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Fig. 2. — Comparison of the TiO temperatures and those from 
the flux integration method (FIM). The different coloured symbols 
represent the FIM temperature when different values of reddening 
were used. The grey points show the FIM fit with the lowest x 
value (see text for details). 

els with the temperatures from the TiO fits but which 
have had the extinction tuned to give the best fit to the 
whole SED. Clearly, the fit to the NIR spectrum is still 
poor, due to the excess of flux around the H-hump. 

In the case of the SED fits, models with higher values 
of Teff are required to fit the NIR continuum, in which 
the contrast of the H-hump is reduced. This then in turn 
requires higher values of Ay ■ While the fits to the IJHK 
region are good, the BVR region is also well matched 
in most cases, apart from the spectral range where the 
TiO features dominate (~600-800nm). This is despite 
the fact that no part of the spectra blueward of 900nm 
was used to tune the SED fits. Only in one object is 
there a substantial discrepancy between the SED fit and 
the blue flux (LMC 137818). ' 

To check whether the H-hump issue is an artifact of 
poor flux calibration, we compared the 2MASS photom- 



etry of each object with our synthetic photometry com- 
puted by convolving the XSHOOTER spectra with the 
2MASS filter profiles. No object has any systematic off- 
set for the i/-band flux compared with J and Ks, and 
the two sets of photometry agree to within the errors. 
Some objects do have 3ct differences in all NIR bands, 
though these differences are not correlated with any stel- 
lar property. We attribute these differences to photomet- 
ric variablity of some of the stars in our sample. 

The differences in the Toff values from the two fitting 
methods are illustrated in Fig.[TJ A systematic difference 
in temperatures is clearly seen. In the case of the TiO 
fits, a systematic offset is seen between the average tem- 
peratures of the stars in the two galaxies. For the SED 
fits, the average Tcff values are roughly the same for both 
galaxies within the errors, 4170±170K and 4030±90K 
for the LMC and SMC respectively. The offset of the 
LMC and SMC stars in Fig. [1] is due entirely to the 
metallicity dependence of the TiO absorption strengths. 

We note that the differences in the Toff and Ay values 
from the two fits cannot be reconciled by changing the 
model metallicities, gravities or microturbulence param- 
eters within sensible boundaries. It is possible to obtain 
similar results for the two analysis methods only if very 
high metallicities are used (i.e. [Z]>0.0 for both LMC 
and SMC), but there is little motivation for such high 
abundances in the Magellanic Clouds (e.g. lTrundle et all 
I2007t ). This wih be discussed further in Sect. O 

As well as the SED Tes values being warmer than the 
TiO fits, the Ay values are also higher (see Table It 
is reasonable to ask whether these higher extinctions are 
believable, since unrealistically high values of Ay would 
put the SED Toff values in doubt. In fact, these val- 
ues of Ay are consistent with that of neighbouring stars, 
plus an extra few tenths of a magnitude for circumstellar 
material, consistent observations of Galactic objects (see 
Sect. l5.11 for further discussion). 

4.2. Results from the flux integration method 

In this method, the extinction Ay is unconstrained. 
Applying a reasonable minimum average reddening value 
towards the Magellanic C louds of E(B-V) = 0.1 mag 
(jZaritskv et al.ll2002l [200l leads to a first resuff which 
is displayed as the red circles in Fig. [5J It appears that 
the temperatures obtained with the FIM-method for this 
minimal value of reddening are on average 200K higher 
than the TiO v alues. Note that we u se the reddening laws 
determined bv I Gordon et al.l (|2003f) . which relate extinc- 
tion Ay and reddening through Ay = RyE{B — V) with 
Ry equal to 2.7 and 3.4 in the LMC and SMC, respec- 
tively. 

Of course, as already indicated by the results obtained 
with the SED-fitting method it is very likely the red- 
dening and extinction towards the sites of the RSGs 
are much higher than the minimal reddening. To ac- 
count for possible reddening values different from O.lmag 
we assume a range of reddening values between E{B — 
y)=O.Omag and O.Smag and determine Toff in the same 
way as above at each reddening (orange, yellow, green 
and blue symbols in Fig. [5]). In addition, we compare 
the dereddened observed fiuxes with the model fluxes 
at three wavelength bands which are free of molecular 
absorption (1.605-1. 702/im, 1.212-1.278Aim, and 1.000- 
1.080/im) and calculate values for each E{B — V). 
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Fig. 3. — Comparison of the MARCS synthetic spectra in the regions of the TiO bands (left) and the H-hump around 1.5-1.6^m (right). 
All models have T(,fl=3700K, logg=0.0, and 5=2km s~^. The metallicities of each model are displayed in the figure legend. 



The minimum of ^ is then used to determine the best- 
fitting E(B — V) and Tcff (grey points in Fig. [5]). We note 
that this extended method aiming at a simultaneous de- 
termination of Teff and E{B — also depends on the 
use of model atmospheres, and is very similar to the SED 
method. As such, these best-fits cannot be considered to 
be entirely independent of the results from the SED fits. 

What these results do show is that, even for minimal 
reddening, the FIM temperatures are always higher than 
the TiO temperatures by 100-200K. Increasing the red- 
dening only serves to increase the discrepancy further. 
We again state that deviations from the standard in- 
terstellar extinction law cannot explain this discrepancy, 
since this can only account for temperature differences 
of ±50K even for the most extreme values of Ry. This 
served to cast further doubt on the reliability of the TiO 
temperatures. 
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5. DISCUSSION 
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5.1. 



Which temperature scale is correct? 

Since the different fitting methods give two different 
temperatures for the stars observed, it is natural to ask 
which, if any, is correct. One major reason to be suspi- 
cious of the TiO temperatures is that they consistently 
overpredict the flux at wavelengths >900nm, dramati- 
cally so in some cases (e.g. LMC 064648, LMC 137818, 
see Fig. [5]) . Though the SED fits overpredict the optical 
flux, this may be plausibly explained by TiO absorption 
in the MARCS models (see later). Conversely, there is 
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Fig. 4. — The amount of reddening deduced from each of the 
TiO band and SED fits, versus the equivalent width of the DIB at 
8620 A. The dotted line sh ows the expected relation in the Galaxy, 
from lMunari et al.l li2008i) 



no obvious physical mechanism that would cause us to 
greatly overestimate the NIR fiux whilst correctly match- 
ing that in the optical. One cannot reduce the extinction 
to rebalance the optical and NIR fiuxes as the Ay values 
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from the TiO fits are already very low. We note that this 
discrepancy was already flagged in L06, who found that 
the V — K colours of RSGs did not match the predictions 
of MARCS models with the same Tcff . 

One other strong piece of evidence that the TiO tem- 
peratures are too low comes from the FIM results. This 
method, which uses the observed flux integrated over all 
wavelengths, is insensitive to model dependent features 
such as the predicted strengths of molecular absorption. 
Instead, it assumes only that the models correctly pre- 
dict the ratios of the total flux to that in a specified 
wavelength interval, which we chose to be a small sec- 
tion of the ii'-band. These results showed that even for 
very low values of reddening the FIM temperatures were 
systematically higher than the TiO temperatures. 

Another line of investigation is to examine the Ay im- 
plied by each fitting method. Here, we consider only the 
TiO and SED results, since the FIM results are not en- 
tirely independent of the latter. For a given fit, to some 
extent there will always be a degree of degeneracy be- 
tween Teff and Ay, since the flux peak can be shifted 
for example to the blue by either increasing Tcff or de- 
creasing Ay'^- As the two methods (TiO and SED) place 
emphasis on different spectral properties, each gives a 
different Ay. Therefore, by finding some other diagnos- 
tic of Ay, we may assess which method is providing the 
most reliable results. 

An independent measurement of extinction may be ob- 
tained from the strengths of the diffuse interstellar bands 
(DIBs). Many DIBs are known to strongly correlate with 
reddening (e.g. ICox et a l. 2006, 2007); however the ma- 
jority of these features are located in the V and R bands, 
meaning that in RSGs they are blended with many in- 
trinsic features. One DIB located in a relatively clean 
spectral region is that at 8620A, identified in the spectra 
of s tars in the RAd ial V elocity Experiment (RAVE) sam- 
ple ([Munari et al.ll2008[ ). In Fig. H we plot the strength of 
this line from each star in our sample against the redden- 
ing E{B — V) obtained from each of the fitting methods. 
We find that the reddenings from the TiO fits appear to 
be poorly correlated with the DIB strengths. The formal 
Pearson coefficient is small and negative (-0.35), imply- 
ing that if anything the data are anticorrelated. There is 
a greater degree of correlation however between the DIB 
strengths and the reddenings from the SED fits. This 
time the Pearson coefficient is positive, though still quite 
low (-1-0.44). A linear fit to the data using the GSFC 
IDL routine fitexy, which considers errors in both vari- 
ables'^, yields a slope and offset consistent with the fit 
to the RAVE stars within the errors (the latter shown 
as a dashed line in Fig. |4|). While this evidence is far 
from conclusive, it does provide circumstantial evidence 
in favour of the SED Tcff values over those from the TiO 
fits. 

We can also estimate the extinction towards each star 
from that of the surrounding stars. Though RSGs may 

^ This applies mainly to the SED method, though it is also 
relevant to the TiO method for earlier spectral types, i.e. when 
the TiO bands are weak. 

^ The routine assumes the errors are symmetric, whereas the 
errors in our E(B-V) values tend to be asymmetric. However, since 
the level of asymmetry is small, and as we are not interested in a 
precise formal fit to the data, we used the larger of the two errors 
in this analysis. 
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Fig. 5. — The circumstellar extinction around each star, found 
by subtracting the local from the total extinction, as a function of 
luminosity. The results of both the SED and TiO fits are shown. 
Note the asymmetric and correlated errors on each datapoint. 

have some extra circumstellar extinction caused by their 
mass-loss history, the surrounding stars should give us 
an idea as to the foreground extinction, both Galactic 
and internal to the star's host galaxy, and so give an 
indication as the the lower limit of the Ay towards each 
object. 

The foreground extinction was measured by looking 
at that towards stars wi thin 1' of our targets from 
IZaritskv et all (|2002l I2004D . We experimented with sev- 
eral methods of interpolating the extinction at the loca- 
tion of our targets, such as binning in right-ascension and 
declination or in radius from the target star. In the end 
we found that simply median averaging the stars within 
0.5' gave results similar to any other method tried, and 
was stable to within ±0.2mags in Ay. 

In Fig. [5] we plot the circumstellar extinction (the to- 
tal extinction from the spectral fits minus the foreground 
extinction from the neighbouring stars) as a function of 
stellar luminosity L (see Sect. 15.41 for details on how we 
measure L). The first thing to note is that the circum- 
stellar extinction implied by the TiO fits is often nega- 
tive - that is, the localized extinction is greater than the 
total extinction measured from the TiO region. In con- 
trast, the total extinction as measured from the SED 
fits is always at least equal to that of the local ex- 
tinction. Some objects, especially the more luminous 
ones, show evidence of substantial circumstellar extinc- 
tion of up to Ay=lmag, consi stent with mid-IR obser- 
vations of Galactic RSGs (e.g. Ide Wit et aLll2008[) . and 
expe ctations calculated from sim plified mass-loss histo- 
ries (jWalmswell fc Eldridgel[20T2l ). In addition, there is a 
separation between the LMC and SMC stars, suggesting 
a possible metallicity dependence on RSG wind density. 

To summarise, we conclude that the Teg measurements 
from the strengths of the TiO lines are unreliable, since 
model fits to these lines overpredict the NIR flux, and 
extinctions implied are anomalously low. The SED tem- 
peratures on the other hand provide good fits to all parts 
of the spectrum apart from the TiO bands. They also 
imply extinctions that are consistent with their neigh- 
bouring stars with a contribution from a circumstellar 
component. In the following section we discuss how the 
mismatch between the SED and TiO results may arise, 
and discuss possible solutions. 
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5.2. Reconcilling the two T^g measurements 

Given this clear discrepancy between the different in- 
dependent temperature measurements, we now explore 
various explanations for this result. We first note that 
the differences in the T^s and Ay values from the two fits 
cannot be reconciled by changing the model gravities or 
microturbulence parameters within sensible boundaries. 

The discrepancy between the two methods may dis- 
appear at higher metallicities. In Fig. |3] we show the 
MARCS model spectra at fixed T^s, logg and ^, for 
SMC-like, LMC-like, Solar, and Solar -HO. 5dex metalhci- 
ties. The strengths of the TiO bands are strongly cor- 
related with [Z], as one would expect. In addition, the 
strength of the H-hump is inversely correlated with [Z], 
as atomic and molecular absorption eat away at the sub- 
peak caused by a minimum in the H~ opacity. These two 
effects may therefore act to bring the two T^s measure- 
ments together: increased TiO absorption would require 
models with higher Toff values to provide a match to TiO 
bands of a given strength; while a reduced H-hump would 
mean that good SED fits would be obtained at lower Tcff 
values. Therefore, potentially one could obtain similar 
results for the two analysis methods if very high metal- 
hcities are used (i.e. [Z]>0.0 for both LMC and SMC). 
However there is little motivation for such high abun- 
dances in the Magellanic Clouds. Young-age metallicity 
tracers such as cepheids, F supergiants and B stars con- 
sistently show average abundances of [Z]lmc = -0-4 and 
[Z]lmc — -0.6, with errors on the means of ±0.05dex, 
and t ypical intrasample disper si ons of less than ifc0 . 02dex 
Ce.g. lAndrievskv et all 120011: iRomaniello et all 120051: 
IKeller fc Woodll2006l: lTru~ et al. 

5.2.1. Adjusting CNO abundances 

Rather than changing the adundances of all metals in 
a uniform way,one might suggest that adjusting only the 
CNO mixture to reflect the products of nuclear burning 
would alter the appearance of the spectrum. The rel- 
ative abundances of TiO, CN, and CO should depend 
on the CNO mix, and all are major sources of opacity. 
Altering the CNO mixture may therefore alter the tem- 
perature structure of the star, changing the strengths of 
the diagnostic absorption features. 

To investigate this, we computed a model with 
Tefr=3700K, [Z]=-0.5, logg^O.O, and ^=2km s^^, with a 
CNO mixture altered to reflect stellar evolution. Specif- 
ically, we reduced C by 0.4dex, and correspondingly in- 
creased N by -|-0.54dex. This is similar to the s urface 
abundances of Galactic RSGs (iDavies etall [20091) . We 
then compared this to a model with the same parame- 
ters other than the CNO mix set to Solar. The evolved 
CNO model had slightly weaker CO and CN features, 
due to the reduced C abundance; while TiO absorption 
was slightly increased due to less O being tied up in CO. 
However, the differences between the two were minimal. 
Analysing the evolved CNO model in the same way as 
the XSHOOTER spectra in this paper, we found that 
this led to Tcff discrepancies of less than 20K. We there- 
fore conclude that deviations from a Solar CNO mixture 
cannot explain the results presented here. 

5.2.2. Convection, granulation, and temperature structure 

As the strengths of the TiO bands are very tempera- 
ture sensitive, inadequacies in the temperature structure 



of the models may lead to a discrepancy in the absorption 
strength as a function of Tgff , giving rise to the effects de- 
scribed in this paper. Indeed, this was suggested in L06 
as an explanation for the V — K discrepancy. Here we 
discuss possible causes of such an effect. 

Convection is known to occur in the envelopes of RSGs, 
and is the physical justification for the inclusion of both 
micro- and macro-turbulence in spectral synthesis mod- 
els. In order to self-consistently model these effects 
rather than including them in an ad-hoc fashion, one 
needs to switch from 1-D hydrostatic to 3-D hydrody- 
na mic stellar models . Such work is still in its early stages, 
but lChiavassa et all (|2011h present a 3-D non-gray model 
of a RSG, along with a synthesized spectrum, and com- 
parisons to 1-D models. These authors show that the 
temperature structure for the 3-D models is subtley dif- 
ferent from that of the 1-D model, resulting in the forma- 
tion zone of the TiO lines shifting outwards to larger radii 
and lower temperatures. In terms of the time-averaged 
SED of the 3-D model, they found that their non-grey 
model had similar TiO band strength to the 1-D model 
with a low Tefr=3400K, but that the SED >l^m looked 
more like the 1-D model with reff=3700K ^. 

In Fig. [B] we show this quantitatively, by fitting the 
Chiavassa et al. <3-D> spectrum with ID MARCS mod- 
els just as in our analysis of the XSHOOTER spec- 
tra of the LMC/SMC stars. We used the same grav- 
ity and metallicity as were computed for the 3-D model 
(log (7=-0.34, [Z]=0.0)^. As with our observations, we 
find that the TiO fit has too much fiux in the NIR, 
while the fitted temperature is lower than the input value 
(Txio = 3580 ± 50K). The SED fit again provides an 
excellent match to all regions 500-2500nm, other than 
in the regions of the TiO bands, and has a tempera- 
ture closer to the integrated fiux of the <3-D> model 
(TsED — 3740±50i4r). Both fitted models are discrepant 
at short wavelengths (<500nm). 

Though a detailed analysis of 3-D model spectra is be- 
yond the scope of this current work, it seems as though 
3-D effects may serve to reconcile the discrepancy in the 
Toff as measured from different regions of the spectrum 
we have highlighted here. In particular, the tempera- 
ture in the outer layers of the 3-D model where the TiO 
lines form is somewhat lower, leading to increased TiO 
absorption. However, the temperature at optical depth 
r w 1, where the continuum forms, is roughly the same 
as in the 1-D model (see their Fig. 5, top-right panel). 
For this reason we conclude that the temperatures mea- 
sured from the continuum regions are more reliable than 
those measured from the TiO bands. 

Finally, we suggest that another potential complicat- 
ing factor could be the difficulty in delineating the stel- 
lar photosphere, the molecular formation zone, and the 
stellar wind. If molecules form in the wind, or close to 
it, then their strengths may depend on wind density. If 

■* IChiavassa et al] 1120111 ) quote 3430K as the temperature of 
their model. However, this is the result of a complex non-standard 
averaging procedure. The integrated flux of their model corre- 
sponds to T(,ft=3700K, which is also that model's temperature at 
r = l. 

^ In their paper, IChiavassa et al] II2011I ) state that due to the 
effect of convective pressure the effective gravity is logg=-0.65. 
However, we found that adjusting the gravity had little effect on 
our fits, other than a change in the opacity below 500nm. 
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Fig. 6. — Analysis of the 3D CO^bold RSG spectrum of IChiavassa et al.l l|2011l ') with the ID MARCS models. The left panel shows 
the 3D spectrum (black) along with the fit to the TiO bands (green), yielding a best-fit T;,ff=3600K. The right panel shows that this 
fit overpredicts the flux at the H-hump, and in the /f-band. Meanwhile, a flt to the NIR continuum regions (magenta) gives a higher 
temperature of 3800K, while underpredicting the strengths of the TiO bands. This mimics the behaviour seen in our spectra of RSGs in 
the Magellanic Clouds (see Fig. |9]l. 

well-known Galactic RSGs, such as fi Cep and a Ori, 
are are too bright to observe with VLT/XSHOOTER, 
which currently is the only instrument capable of simul- 
taneous optical/near- IR spectrophotometry. Meanwhile, 
more distant RSGs suc h as those in the Scutum RSG 
clust ers at 6.6kpc (e.g. iFiger et al] 120061 : iDavies et al.l 
l2007f ) have large foreground extinction of Ay ^15 mag, 
making optical spectrophotometry of the TiO bands ex- 
tremely challenging. 

We have obtained NIR spectroscopy of stars in 
Per OBI with IRFT-|-SpeX (Gazak et al, in prep), and 
combined these data with the optical spectroscopy of 
LOS. We then analysed these stars in the same way as 
our XSHOOTER data. We found very similar results, 
specifically that the SED temperatures were consistently 
4100±150K, whereas the TiO temperatures were much 
cooler and varied with spectral type (agreeing well with 
the results of L05). 

However, there are potential causes for concern with 
these data. Firstly they are non-contemporaneous, the 
optical and NIR spectroscopy taken several years apart. 
Secondly, the SpeX slit was <1" wide, leading to in- 
evitable slit losses and compromising the accuracy of 
the spectrophotometry. Finally, the absolute photomet- 
ric calibration was done with 2MASS, in which the stars 
are saturated and have photometric errors of ±0.2mags. 
Since the data quality does not match that of our excel- 
lent XSHOOTER data, we have chosen not to present 
the results of that analysis here. 

5.4. Implications for stellar evolution 

To see how the difference in Toff affects the stars' posi- 
tions in the H-R diagram, we determine the luminosities 
of the stars for the TiO and SED Toff estimates. Here, 
we use distance moduli of 18.47 and 18.95 for the LMC 
and SMC respectively (averages of a number of measure- 
ments, taken from the NASA Extragalactic Database). 

For the TiO temperatures we first take the F-band 
magnitude of each star from our synthetic photometry. 
Wc then correct for the Ay determined from the TiO 
fits, and apply the bolometric corrections from L06. 
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Fig. 7. — HR-diagram of the stars in our sample. The magenta 
points show the locations of the stars when the TiO temperatures 
and extinctions are used, while the blue points use the values de- 
rived from the SED fits. Note that the errors on each quantity 
are correlated and asymmetric. The grey lines show the Geneva 
mass tracks for stars with initial rotation 400 km s~^ and LMC-like 
metallicity. 

the ma ss-loss rates of RSGs are a function of their lumi- 
nosity (jBonanos et al.|[2010f ). then the strengths of the 
molecular absorption bands would be highly sensitive to 
luminosity as well as temperature. If mass-los s rates are 
also [Z] dependent (jMauron fc Josselinl[20nl see also 
Fig. E] in this paper), then this may serve to increase this 
sensitivity of the band strengths to luminosity. 

5.3. Comparison with Galactic RSGs 

Since there is a known dependence of the TiO band 
strengths on metallicity, it is interesting to ask how a 
sample of Galactic RSGs would compare to the LMC and 
SMC stars studied here. However, obtaining data com- 
parable to that presented here is problematic. The most 
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For the SED fits, the luminosity was obtained by inte- 
grating under the SED. The flux in the regions of high 
telluric extinction was interpolated using the best-fitting 
MARCS model for a given star. The flux longward of 
2.5/im was estimated by logarithmic interpolation be- 
tween a star's flux at 2.3/j,m (i.e. before the CO bandhead 
absorption kicks in) and the four Spitzer/IRAC fluxes at 
3.6/im, 4.5^m, 5.8/j,m and 8.0/im. The total flux outside 
our observed region of 400nm-8.0/im was c onsidered neg- 
ligible . The flux was dereddened using the iGordon et al.l 
([2003') extinction law, and assuming that the correction 
for wavelengths greater than 2.5/im was effectively zero. 
The bolometric flux was then calculated by integrating 
under the dereddened SED. 

The H-R diagram for the stars studied here is shown 
in Fig. [71 When the temperatures and luminosities are 
updated to reflect the results of the SED analyses, we find 
that the stars move to warmer temperatures at roughly 
constant luminosity. The reason that the luminosities are 
largly unchanged is that, when analysing the SEDs, we 
find that both Toff and Ay are higher than the results 
of the TiO fits. This means that, in the TiO fits, the 
smaller extinction correction is cancelled out by a larger 
bolometric correction. 

In terms of the location of the stars in the H-R diagram 
in relation to the predictions of stellar evolutionary the- 
ory, the SED temperatures bring the RSGs into much 
better agreement with the Geneva rotating model calcu- 
lations. It was pointed out by LOG that the TiO temper- 
atures were slightly cooler than the coolest point (i.e. the 
Hayashi limit) of the Geneva stellar tracks. We note how- 
ever that the other contemporary stellar e volution codes, 
such as the Bonn or C ambridge codes ([Eldridge et al.l 
120081 iBrott et al.l 1201 ll ) predict that the Hayashi limit 
at LMC metallicity should be cooler than that of the 
Geneva models, at around 3700K. It is therefore worth 
discussing the nature of the Hayashi limit of RSGs a little 
further. 

The Hayashi limit is the maximum size to which a 
star can grow whilst remaining in hydrostatic equilib- 
rium, and depends on a number of factors, but is mainly 
governed by the opacity of the envelope (i.e. metallic- 
ity) and the conv ective mixing length parameter a (e.g. 
iLucv et al1ll986l ). In most stellar codes, a is tuned to 
reproduce the observed properties of the Sun, and is as- 
sumed to be the same for stars of all masses and metallic- 
ities. Some authors have argued that a single value of a 
cannot reprod uce the colour-magn itude diagrams of glob- 
ular clusters ([Ferraro et al.l l2006f) , while there is some 
deb ate about whether a may be depende nt on metallic- 
itv (ISalaris fc Cassisil [1996: Keller 1999: iPalmieri et al.l 

mm- [ 

Given this discussion, it is perhaps surprising that 
though the Bonn, Cambridge and Geneva codes all use 
the same mixing length, they predict different tempera- 
tures for the Hayashi limit of RSGs. The reasons for this 
are unclear. The solution to this problem will ultimately 
come from reliable temperature measurements of RSGs. 
The results presented here are therefore a significant step 
towards understanding the nature of convection in cool 
stars. 

5.5. Bolometric corrections, luminosity calibrations, 
and application to supernova progenitors 
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Fig. 8. — The pre-explosion photometry of SN2008bk's progeni- 
tor, and the best-fit model with T^ff =4200K. The photometry from 
IMattila et a h (2008) is represented by the black triangles, while the 
synthetic photometry from the best-fitting MARCS model is shown 
by the blue circles. The full SED of that model is shown by the 
cyan line. 



TABLE 3 

Mean bolometric corrections and standard deviations for 
RSGs IN each of the two Magellanic Clouds. 



Galaxy 


LMC 


SMC 


BCk 


2.69 ± 0.12 


2.69 ± 0.06 



TABLE 4 

Parameters for use with Eq. Q to derive bolometric 
luminosities from single-band photometry of rsgs. 



Band 


a 


b 


(Tlog(L) 


V 


3.12 ± 0.06 


-0.29 ± 0.01 


0.12 


R 


2.44 ± 0.07 


-0.34 ± 0.01 


0.09 


I 


1.90 ± 0.08 


-0.37 ± 0.01 


0.06 


J 


1.30 ± 0.09 


-0.39 ± 0.01 


0.03 


H 


0.97 ± 0.10 


-0.40 ± 0.01 


0.04 


K 


0.90 ± 0.11 


-0.40 ± 0.01 


0.04 



5.5.1. Bolometric corrections 

The uncertainties that we highlight in the RSG tem- 
perature scale are extremely important for the study 
of SN progenitors. Object brightnesses and colours in 
pre-explosion images are used to place the progenitor on 
the HR diagram, and by comparison with stellar mod- 
els determine the progenitor mass. Clearly, such work 
is sensitive to both the progenitor object's extinction 
and, through the bolometric correction, the calibration 
of colour with effective temperature. 

We calculate the BC for each star in our sample from 
the bolometric flux, as determined from the SED fits in 
Sect. 15.41 the flux from the synthetic photometry (Sect. 
12. 2p once corrected for extinction (from the SED fits), 
and assuming an absolute bolometric magnitude for the 
Sun of 4.77. 

We find that, as predicted, the BCs at various bands 
are correlated with temperature. However, we have 
shown here that the temperatures of RSGs, as deter- 
mined from the continuum regions of the SED, are uni- 
form to within iblOOK. This means that the amplitude 
of variations in BC are small, particularly in the if-band 
where we are close to the flux peak and where there is lit- 
tle contribution from molecular absorption. Hence, to a 
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good approximation we can assume that for RSGs BC^ 
is constant to within '^0.1 mag. 

In Table [3] we hst the average BCs for the RSGs in 
both the LMC and SMC in the Ks band. The BCs for 
both galaxies are consistent with one another to within 
the errors. This has the important implication that we 
can determine the bolometric apparent magnitude of a 
RSG - and if the distance is known, the luminosity - 
from only one photometric data-point. 

As is to be expected, the BC^ values presented here 
are much changed compared to those quoted in L06. This 
is caused by the incorrect assumption that the models 
which best-fit the TiO bands will accurately predict the 
flux at all wave lengths. In fact, t he LMC BCk scale is 
closer to that of lElias et al.l (|1985[ ). which relied on aver- 
aged instrinsic colours for M supergiants and trapezoidal 
integration of optical/near-IR broad-band photometry. 
Our results for the SMC however are greatly different 
to those of Elias et al., since those authors did not have 
accurate intrinsic colours for SMC RSGs. 

5.5.2. Luminosity calibrations 

To be able to determine the luminosity of the star from 
a given photometric point, we must first measure and 
subtract the extinction. We separate this extinction into 
two components, interstellar (IS) and circumstellar (CS). 
The IS extinction is measured from the surrounding stars 
(Sect. [5T|) . whereas the total (IS-I-CS) is measured from 
the SED fits. Since the CS component is typically small 
(Av^l)^ in the infrared bands this correction will be 
minor. Therefore, if the IS component can be estimated, 
one can estimate the luminosity of a RSG from i^-band 
photometry alone. Furthermore, since the CS extinction 
is somewhat correlated with luminosity (see Sect. 15. 1|) . 
it is possible to calibrate this effect out and estimate the 
luminosity from any pass-band. 

To do this, we plotted the magnitudes of each star 
though a given pass-band, corrected for foreground ex- 
tinction, against the luminosities from the SED fits. We 
then used the IDL procedure linfit to determine the 
coefficients of the linear fits a and b to these relations, 
which we list in Table |4l The bolometric luminosity of a 
RSG can then be determined from, 

log{L/LQ) = a + b{mx - (2) 

where m\ is the apparent magnitude at a given pass- 
band after correction for interstellar extinction, and /x is 
the distance modulus. In the cases of IJHK, b is close to 
0.4(=2.5~^), which follows simply from the magnitude 
definition. At V and R bands, b begins to diverge from 
this value, which is due to the fiux at these bands be- 
ing sensitive to CS extinction and TiO absorption, both 
of which are luminosity dependent. The final column 
of Table H] shows the r.m.s. standard deviations of the 
residuals between the luminosity from the SED fits and 
those determined from Eq. ([2]). In practice, the dominant 
sources of error on log(L/LQ are the uncertainties on a 
and b, which corresponds to ~0.1dex at K and ~0.20dex 
at V. 

We note that we found no appreciable differences in 
the results when the sample was divided into stars from 
each of the LMC and SMC. We therefore combined all 
the stars to give greater signal-to-noise. In future, using 



larger samples of stars, these relations could be improved 
upon by analysing the stars from the two galaxies seper- 
ately to take account of metallicity effects, increasing the 
accuracy. 

5.5.3. Application to pre-explosion photometry of Type II-P 

supernovae 

To assess the potential impact on the field of SN pro- 
genitors we re-analyse the pre-explosion photometry of 
the Type II-P supernova 2008bk using our SED temper- 
atures. This object was chosen as it has LMC-like metal- 
licity, consistent with the objects studied in this cur- 
rent work, and pre- explosion photometry in four bands 
(jMattila et al.l 120081 ). rather than the usu al one or two 
which is the case for most progenitors (jSmartt et al.l 
I2OO9I) . In their paper, Mattila et al. fit the photometry as 
an M4 super giant, Av=l.Q± .5, an d using the tempera- 
ture scale of iLevesaue et al.l (|2005[) . BCk = 2.9 ± 0.1. 
This led to a bolometric luminosity of log{L/LQ) = 
4.6 ±0.1, and by comparison to stellar models, an initial 
mass for the progenitor star of Q^^M©. They note that 
the SED may also be fit by a MO supergiant and a higher 
extinction of Av=3.0. 

To re-analyse the pre-explosion photometry of 
SN2008bk, we restrict the model Tcs to between 
4130±150K, consistent with our results for the LMC, 
and use an LMC-like metallicity of [Z]=-0.4. We again 
use log5=0.0 and ^ = 2.0km s~^, noting that these lat- 
ter two parameters make little difference to the results. 
For each temperature we determine the Ay required to 
produce the best fit. 

We find that the photometry can be fit very well by 
models in this temperature range, requiring an extinc- 
ti on of Av=2A±0A. Following the analysis procedure 
of Ma ttila et al.l ()2008D . we use the object's i^-band pho- 
tometry (if =18.34±0.07), the bolometric correction ap- 
propriate for the LMC (BCk = 2.69 ±0.12), the extinc- 
tion at K (=At/x0.11 iRieke fc Lebofskvl [l985l ). and a 
distance modulus of 27.9±0.2. This gives us a bolometric 
luminosity for the progenito r of log(L/£(^) = 4. 75 ±0.10, 
slightly higher than the Mattila et alJ (|2008[ ) estimate. 

We can also estimate Lboi using our emprical relations 
in the previous section. We use the object's K-hand 
photometry, since this is the least sensitive to extinc- 
tion. Assuming negligible foreground extinction, we find 
log(i/L0) = 4.72 ±0.14. Mattila et al. suggest that the 
foreground Galactic extinction to the host galaxy is low 
(Ay~0.06), but that a comparison to a SN with a similar 
light-curve suggests that there may be around O.Smag of 
extinction at V internal to that galaxy. This is a negli- 
gible amount of extinction at K, and would increase the 
luminosity of the progenitor by only ^^O.Oldex. 

To convert this luminosity into an initial stellar mass, 
we compare to the stellar evolutionary tracks at LMC 
metallicity. We use the tracks of both the Bonn 
(Bro tt et aLll2011[ ) and Geneva groups. Our assumed 
temperature for the progenitor star means that its lo- 
cation in the HR-diagram does not coincide with the 
terminal point of any of these tracks, which evolve to 
temperatures of ~3500K. However, as noted earlier, this 
terminal T^g is set by the adopted value for the mixing 
length in the stellar structure models, which is poorly 
constrained for massive stars. In order to estimate the 
initial mass of the progenitor, we follow the methodology 
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of lSmartt et al.l (|2009l . and references therein) and com- 
pare the star's luminosity to the terminal luminosities of 
the mass tracks, which corresponds to the luminosity at 
the end of He burning. In this way, we find a progeni- 
tor mass for SN2008bk of M = U^IMq using the Bonn 
tracks, and M = 11 ± IMq with Geneva t racks. This is 
^30% higher than the previous estimate of lMattila et al.l 
(|2008l ). though within their errors. 

We note tha t our extincti o n is slightly higher than 
that derived bv lMattila et al.l (|2008), and that those au- 
thors argued against a higher Ay solution due to the lack 
of signs of circumstellar extinction in the SN spectrum. 
However, it is entirely plausible that for a relatively small 
amount of circumstellar extinction close to the star, this 
material would be dest royed al most instantly by the SN 
explosion. Indeed, Fraser et al.l ([2012. ) make a similar ar- 
gument for the progenitor of SN2012aw, for which multi- 
band pre-explosion photometry was also available. 

In summary, we find that if we take the temperatures 
of RSGs as derived from fits to the near-IR SED, the 
extinction increases, while the magnitude of the bolo- 
metric correction decreases. Though these two effects 
work in opposite directions, in the case of SN2008bk the 
net effect is for the progenitor star's inferred bolometric 
luminosity to increase by 0.14dex. Therefore, the ini- 
tial mass of the p rogenitor is higher compare d to that 
derived using the iLevesaue et al l (|2005l . I2006D temper- 
ature scale. If the change in temperature scale affects 
the properties of progenitors in the same way, this may 
in part exp lain the so-called 'Red Supergiant problem' 
(jSmartt et al. 2009), whereby the progenitor masses of 
Type II-P SNe tend to be lower than those expected 
from initial mass function arguments. This is similar to 
the suggestion of iWalmswell fc Eldrided (|2012f ) that the 
'missing' SN progenitors could be RSGs which have suf- 
ficient circumstellar extinction to push them below the 
detection limits of the pre-explosion imaging. 

5.6. Concluding remarks on the temperatures of 
Red Supergiants 

In the course of this work, we have shown that for 
RSGs the definition of a characteristic temperature is a 
non-trivial task. The definition of effective temperature. 
Toff = i/(47ri?^cr), whilst unambiguous for stars such 
as the Sun with well-defined photospheres, can be prob- 
lematic to define for stars with extended atmospheres. 
Simply defining the photosphere to be where the opti- 
cal depth T = 2/3 is insufficient, since this value can de 
reached at very different radii depending on the observed 
wavelength. Any flux-averaged value of r will therefore 
be sensitive radiative transfer processes in the stellar at- 
mosphere. 

Another problem is finding accurate and robust spec- 
tral diagnostics of T^s- For M supergiants, it has always 
been assumed that the strengths of the TiO bands, which 
define the spectral classification sequence of these stars, 
were useful diagnostics of effective temperature. How- 
ever, any diagnostic line is sensitive only to the local 
temperature in the layers where the line forms. In the 
case of the TiO lines, these form high in the atmosphere, 
and so are very sensitive to the radial temperature struc- 
ture. As we have shown in our comparison to 3-D hy- 
drodynamic models, two stars with the same luminosity 
and radius can have different TiO strengths if the tem- 



perature structure is different. However the continuum, 
which forms at much deeper layers, is indistinguishable. 
Though there is potentially a degeneracy with redden- 
ing, this degeneracy can be broken by considering the 
magnitude of the H~ opacity minimum. We therefore 
conclude that the most robust diagnostic of Tcff in RSGs 
is the line-free continuum in the IJHK region, including 
the iJ-hump. 

5.6.1. A connection between spectral type and evolutionary 

stage? 

Finally, it is worth posing the question that if the tem- 
peratures of RSGs are all roughly constant, what is the 
physical interpretation of the spectral type sequence? It 
has long since been assumed that the MK classification 
system was associated with a decreasing temperature 
scale from O to M types, with the subtypes from late-K 
onwards being determined from the strengths of the TiO 
bands. This assumption may well be valid for M giants 
and dwarfs, in which the pressure scale-heights and con- 
vective cells are much smaller with respect to the size 
of the star. However, in the case of RSGs, our results 
from SED fitting show that the temperatures are not 
strongly correlated with the strengths of the TiO bands 
(and therefore spectral type). Inspection of Fig.jTjreveals 
that the Ttio values (a proxy for spectral type) do seem 
to be correlated with luminosity, si milar to the Ca II 
lines used in extragalactic studies fe.g. lHumphrevs et al.l 
Il986( l. Therefore, in the M supergiant domain, it seems 
the spectral types of stars are highly sensitive to lumi- 
nosity as well as temperature. 

In discussing this further, we first make the assumption 
that the TiO lines, which form far above the continuum- 
forming layers, are sensitive to the star's pressure scale- 
height h. This parameter goes as /i ~ T/g, from which 
it is trivial to show that h ^ L/T^, or if normalised 
to the size of the star, h/R ^ L'^-^/T (assuming con- 
stant mass). Even using the previous temperature scale 
of L05, the range in temperatures of RSGs is quite nar- 
row, -3600-4400K (or ±10%). By contrast, evolution- 
ary models show that one expects a RSG to increase 
in L by a factor of ^2 as it evolves (see e.g. Fig. 8 in 
iMevnet fc Maedeil l2000f ) . Therefore, as a star ascends 
the RSG branch, its scale-height will increase, driven pri- 
marily by the increasing luminosity. 

From this we suggest that RSGs with the earliest spec- 
tral types may be those objects which have just arrived 
in the RSG phase, having evolved from the blue. The 
time-averaged spectral type of a RSG will then move 
steadily to later types as the star gradually increases in 
luminosity as it evolves. This is consistent with observa- 
tions of coeval clusters of RSGs, where all RSGs should 
have approximately the same initial mass, in which a 
correlatio n is found between MK sp ectral type and mid- 
IR excess ijCohen fc Gaustadlll973D . Larger excess emis- 
sion is indicative of larger amounts of circumstellar ma- 
terial, which suggests higher mass-loss rates (associated 
with higher luminosity, and hence more evolved objects), 
and/or a star which has spent longer in the RSG phase. 

5.6.2. Future work 

In the present study, we have highlighted the differ- 
ences and similarities between the synthesised spectra of 
1-D and 3-D models. Concentrating on these similarities. 
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we have redefined the temperatures of RSGs using the 1- 
D models, and provided results which should be robust 
to 3-D effects. In the future we intend to make a more 
thorough analysis of this by repeating this study with a 
grid of 3-D models. The switch from 1-D hydrostatic to 
3-D hydrodynamic models is non-trivial as it requires a 
dramatic increase in computational effort, nevertheless 
work is currently underway to construct such a grid. 

With a suite of 3-D models we will also address the 
effect of any departures from spherical symmetry. Should 
the surface flux of the star be non-uniform, e.g. due to 
hot-spots and cool-spots, the relation of Eq. ^ would 
no longer hold. This would then cause the observer to 
underestimate the size of the star for a given temperature 
and bolometric flux. This is important with respect to 
comparisons with evolutionary models, which compute 
an effective temperature from the size and luminosity of 
a star assuming spherical symmetry. 

6. SUMMARY & CONCLUSIONS 

We have re-appraised the temperatures of Red Su- 
pergiants (RSGs) in the Magellanic Clouds using 
XSHOOTER flux-calibrated spectroscopy from 400- 
2500nm. We simultaneously flt effective temperature Tcff 
and extinction Ay by fitting the spectra with a grid of 
1-D spherical MARCS model atmospheres. 

When fitting the VRI spectral region, which is domi- 
nated by the TiO absorption which has historically been 
used to define these stars' temperatures, we find temper- 
atures much lower than from fitting the continuum SED 
in the IJHK region. After exploring many avenues of ev- 
idence, we conclude that the SED temperatures are the 
more reliable, for the following reasons: 

• Models which are tuned to match the TiO bands 
overpredict the flux longward of ^900nm, dramat- 
ically so in some cases. The SED flts on the other 
hand provide good to all regions of the spectrum, 
apart from the regions where the TiO bands dom- 
inate. 

• The TiO temperatures Txio a-re lower than those 
measured from comparing the integrated fluxes of 
the data with those of models (the flux integration 
method, FIM), which is the least model-dependent 
method of determining temperature. Even if we 
assume minimal reddening, the FIM temperatures 
TpiM are higher than those from the TiO fits. The 
FIM temperatures do agree with the SED fits if we 
assume a modest amount of extinction, an entirely 
reasonable assumption. 

• The implied total line-of-sight extinctions from the 
TiO fits are in most cases lower than those of 
nearby stars, even though RSGs are expected to 
have higher extinction due to circumstellar mate- 
rial. By contrast, the SED fits imply extinctions 
which are in agreement with or in excess of that of 
the nearby stars. 

• The extinction from the TiO fits seems to be com- 
pletely uncorrelated with the strengths of the dif- 
fuse interstellar band at 862nm, whereas the extinc- 
tion from the SED fits appears to show a degree of 
correlation. 



To reconcile these two measurements of temperature, 
we considered the effects of overall metallicity, CNO 
mixture, and dynamical effects (convection and granu- 
lation). Of these we conclude that the latter is most 
likely cause. Recent studies have shown that 3-D hydro- 
dynamical models of RSGs have a significantly different 
radial temperature profile in the upper layers of the at- 
mosphere where the TiO lines form. The temperature 
structure at deeper layers, where the continuum forms, 
is largely unchanged. For this reason we consider the 
IJHK continuum region and the opacity minimum 
at 1.6/zm to be robust temperature indicators. 

Adopting the SED temperatures, we find a number 
of interesting effects. Firstly, the temperatures of all 
RSGs in both Magellanic Clouds are roughly uniform, 
4150±150K. This brings the observed Hayashi limit for 
massive stars into agreement with those of Geneva mod- 
els, and somewhat warmer than those of Bonn and Cam- 
bridge models. This also means that the strengths of the 
TiO bands, which define the spectral types of M super- 
giants, may also strongly depend on luminosity. From 
this we suggest that the spectral type of a RSG may be 
indicative of its evolutionary state, with more evolved 
RSGs having later spectral types. 

Secondly, the circumstellar extinction of RSGs is 
weakly correlated with luminosity, consistent with recent 
findings that more luminous stars have stronger winds. 
We also see a separation between LMC and SMC stars, 
implying a metallicity dependence on wind strength. 

Thirdly, as there is very little variation in tempera- 
ture of the RSGs, their bolometric corrections are all 
approximately the same. This means that, if we can ac- 
count for extinction, we can determine a RSG's bolomet- 
ric luminosity to a high degree of precision from a single 
photometric point. This has applications for the study 
of supernova progenitors, where pre-explosion imaging 
is often limited to one or two photometric bands. We 
have shown that the circumstellar extinction is very low 
in the near-infrared, so if the foreground extinction is 
known (e.g. from neighbouring stars) the luminosity can 
be measured to a precision of itO.ldcx from only one 
photometric band. In the optical ( VRI) region, circum- 
stellar extinction becomes more significant, as does the 
effect of TiO absorption on the broad-band flux. How- 
ever, since both these effects are luminosity dependent, 
we have been able to provide a calibration which recovers 
the luminosity from optical photometry to a precision of 
±0.2dex. 
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Fig. 9. — MARCS model fits to the XSHOOTER data. The left panels show the observations (black) and the best fitting model (green), 
along with the fit parameters in the top right. The right panels show the same data over a broader wavelength range (black), the best SED 
fit (magenta), the TiO bands fit from the left panel (green), and the best fitting model with the T^f{ from the TiO bands fit but which has 
had the extinction tuned in order to give the best fit to the whole SED (grey). Overplotted in orange are the 2MASS photometry for each 
star. 
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Fig. 9. — Continued. 
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Fig. 9. — Continued. 
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